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ABSTRACT 


An  investigatic -  was  performed  to  evalua-e  the  accuracy 
of  numerical  mode  -  of  near-fields  for  antennas  on  or  near 
surfaces  using  -  Numerical  Electromagnetic  Code  (NEC) . 
Average  power  gain  and  input  impedance  were  calculated,  for 
two  models.  The  first,  a  dipole  antenna  located  inside  a 
cube-shaped  wire-grid  box  of  1  meter  sides  was  evaluated  for 
a  wide  range  of  frequencies  in  free  space.  The  second,  a 
monopole  antenna  mounted  on  the  top  of  two  cube  shaped  boxes 
(a  wire  grid  and  a  surface  patch  box)  of  .1  meter  sides  over 
a  perfectly  conducting  ground  was  evaluated  from  1  to  1.4 
GHz.  The  monopole  was  positioned  at  the  center,  at  an  edge, 
and  at  the  corner  of  the  box  top.  For  the  dipole  in  the 
box,  the  average  gain  and  input  impedance  are  presented  and 
evaluated.  For  the  monopole  on  the  box,  calculated  results 
are  compared  to  measured  admittance  values.  The  NEC  wire- 
grid  model  results  compare  closely  to  the  measurements,  but 
for  surface-patch  models,  only  one  position  of  the  monopole 
yielded  satisfactory  correlation  to  measurements.  Recommen- 
dations for  improved  numerical  modeling  perfonnance  are 
made. 
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I.  INTRODUCTION 

A.  GENERAL  BACKGROUND 

Because  shipboard  operations  are  carried  out  within 
fixed  (small)  distances  from  transmitting  antennas,  the  Navy- 
has  a  unique  and  long-standing  operational  problem:  the 
radiation  from  these  antennas  can  be  dangerous  to  personnel, 
ordnance,  fuel,  and  electronic  equipment.  This  is  called 
RADIATION  HAZARD  (RADHAZ) . 

Accordingly,  the  Navy  has  been  pursuing  the  study  of 
near  fields  of  antennas  for  a  number  of  years.  However,  the 
near  field  structure  is  complex,  and  previous  theoretical 
analysis  has  been  practical  only  for  simple  antennas  in 
uncomplicated  geometrical  settings.  With  the  advent  of  the 
modern  high-speed  computer,  approximate  solution  techniques 
such  as  the  method  of  moments  used  by  the  Numerical  Electro- 
magnetic Code  have  become  practical. [Ref.  1] 

B.  NEED  FOR  STUDY 

With  the  added  emphasis  on  near  field  (due  to  RADHAZ)  by 
EPA  and  USN  demands,  we  must  be  more  accurate  in  our  near 
field  predictions  to  assure  that  personnel  and  equipment 
can  continue  to  function  safely  in  modern  shipboard 
RF  environments. 


C.   SCOPE  AND  LIMITATIONS 

Two  theoretical  models  of  cube  shaped  boxes  (similar  to 
those  found  on  structures  on  shipboard  systems)  are 
examined.  These  boxes  are  modeled  numerically  and  the 
parameters  of  one  of  the  models  is  compared  to  measured 
values  from  Reference  2.  Modeling  techniques  are  varied  to 
assess  their  impact  on  accuracy.  Guidelines  are  developed 
for  maximum  accuracy.  Also,  it  is  known  that  the  Numerical 
Electromagnetic  Code  has  been  found  to  suffer  loss  of 
accuracy  in  VLF  applications  involving  electrically  small 
antennas  [Refs.  3  and  4].  This  may  affect  the  modeling 
efforts  of  this  thesis. 


II.   NEAR  FIELD  THEORY 

A.   EXTERIOR  FIELDS  OF  RADIATING  ANTENNAS 

Figure  2 . 1  describes  the  regions  into  which  the  exterior 
fields  of  a  radiating  antenna  are  commonly  divided.  The 
antenna  radiates  into  free  space  as  a  linear  system  with  the 
single-frequency  time  dependence  of  e^"^*^)  . 

The  FAR-FIELD  region  extends  to  infinity  and  is  that 
region  of  space  where  the  radial  dependence  of  electric  and 
magnetic  fields  varies  approximately  as  e(~^^^)/r.  The 
inner  radius  of  the  far  field  can  be  estimated  from  the 
general  free-space  integral  for  the  vector  potential  and  is 
usually  set  at  A  +{2D^/\)  .  (The  added  \  covers  the  pos- 
sibility of  the  maximum  dimension  D  of  the  antenna  being 
smaller  than  a  wavelength.  In  other  words,  the  RAYLEIGH 
distance  (Z)-'-  should  actually  be  measured  from  the  outer 
boundary  of  the  reactive  near  field  of  the  antenna.)  For 
the  main  beam  direction,  the  RAYLEIGH  distance  can  sometimes 
be  reduced. 

The  free-space  region  from  the  surface  of  the  antenna  to 
the  far  field  is  labeled  as  the  NEAR-FIELD  region.  It  is 
divided  into  two  subregions,  the  reactive  and  the 
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Figure  2.1   Exterior  fields  of  radiating  antennas  [Ref.5] 


radiating  near  fields.  The  reactive  near-field  region  is 
commonly  thought  to  extend  about  A/27r  from  the  surface  of 
the  antenna,  although  experience  with  near-field  measure- 
ments indicates  that  a  distance  of  a  wavelength  (A)  or  so 
would  be  a  more  reasonable  outer  boundary  to  the  reactive 
near  field. 

The  reactive  near  field  can  be  defined  in  terms  of 
planar,  cylindrical,  or  spherical  modes.  Unfortunately,  the 
reactive  fields  of  spherical  (or  cylindrical)  multipoles  are 
not  identical  to  the  plane-wave  evanescent  fields  of 
multipoles.  Thus,  a  simpler,  physically  appealing  and  less 
ambiguous  method  of  defining  the  reactive  region  of  antennas 


relies  directly  on  Poynting's  theorem  and  the  vector 
potential.  It  can  be  shown  that  the  contribution  to  the 
reactive  part  of  the  input  impedance  of  an  antenna  from  the 
fields  outside  a  surface  surrounding  the  antenna  is 
proportional  to  the  imaginary  part  of  the  complex  Poynting 
vector  integrated  over  the  surface.  Thus,  wherever  the 
phase  of  the  electric  and  the  magnetic  field  vectors  are 
near  quadrature,  the  Poynting  vector  will  contribute  mainly 
to  the  reactive  part  of  the  input  impedance.  Taking  the 
curl  of  the  vector  potential  integral  once  to  get  the 
magnetic  field  and  twice  to  get  the  electric  field,  shows 
that  the  phase  of  the  electric  and  magnetic  fields  may  be 
near  quadrature  in  regions  within  approximately  a  wavelength 
(X)  of  the  antenna.  Consequently,  the  region  within  about  a 
wavelength  of  the  physical  antenna  is  referred  to  as  the 
REACTIVE  NEAR  FIELD. 

Beyond  a  distance  of  about  a  wavelength,  the  electric 
and  magnetic  fields  tend  to  propagate  predominantly  in 
phase,  but,  of  course,  do  not  exhibit  e  (~J^^)/r  dependence 
until  the  far  field  is  reached.  This  propagating  region 
between  the  reactive  near  field  and  the  far  field  is  called 
the  RADIATING  NEAR  FIELD. 

Finally,  the  optical  terms  "Fraunhofer  region"  and 
"Fresnel  region"  are  sometimes  used  to  characterize  the 
fields  of  antennas.  The  term  FRAUNHOFER  REGION  can  be  used 
synonymously  with  the  far-field  region  or  to  refer  to  the 


focal  region  of  an  antenna  focused  at  a  finite  distance. 
The  FRESNEL  REGION,  which  extends  from  about  (D/2\)'^/^  D/2 
+  X  t°  "t^®  f^^  field,  is  the  region  up  to  the  far  field  in 
which  a  quadratic  phase  approximation  can  be  used  in  the 
vector  potential  integral.  The  Fresnel  region  is  a  sub- 
region  of  the  radiating  near-field  region. [Ref .  5:pp.  33-34] 

B.   DIFFICULTIES  WITH  NEAR  FIELD  CALCULATIONS 

The  NEAR  FIELD  is  more  difficult  to  calculate  than  the 
far  field.  When  calculating  the  near  zone  radiation,  the 
antenna  no  longer  appears  as  an  infinitesimal  point  (as  is 
the  assumption  for  far  field) ,  but  rather  we  are  conscious 
that  the  energy  coming  toward  us  comes  from  separated 
locations.  Also  the  terms  of  the  field  expressions  with 
powers  of  1/r  (r  is  the  distance  from  the  origin  of  the 
antenna  to  the  field  point)  greater  than  one  are  here  more 
important  than  the  1/r-dependent  terms.  The  complex 
Poynting  vector  ^(E  x  H*)  will  then  contain  terms  with  a 
power  of  (1/r)  in  addition  to  the  radiation  field  term.  It 
results  that  these  terms  are  purely  imaginary,  indicating 
reactive  power  and  energy  oscillating  in  and  out.  The  near 
field  is  restricted  to  charge  and  current  density  and  far 
field  only  to  current.  All  these  difficulties  make  the  near 
field  calculation  more  imprecise. 

In  the  past,  accuracy  of  near-field  calculations  was  not 
a  problem;  there  was  no  practical  application  to  which  the 
energy  could  be  applied.   However,  in  recent  years  there  has 
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been  an  escalating  interest  in  the  close-in  problem, 
especially  in  the  Navy  where,  in  general,  antennas  mounted 
on  a  ship  will  couple  into  metallic  parts  of  the  ship.  This 
means  that  strong  RF  currents  can  be  excited  on  closely 
coupled  conducting  surfaces,  which  in  turn  re-radiate.  The 
total  near  field  is  the  (vector)  sum  of  the  fields  radiated 
by  the  antenna,  the  primary  antenna,  and  the  nearby 
conducting  objects.  Re-radiation  is  enhanced  whenever  the 
secondary  scatterers  are  similar  in  size  to  the  prymary 
antennas  or  whenever  resonant  length  conducting  paths  and 
spacing  occur.  Because  of  this  the  near  field  produces 
RADIATION  HAZARD  (RADHAZ)  that  can  affect  personnel, 
ordnance,  and  fuel  on  ships. 


III.  NUMERICAL  MODELS  OF  SURFACE  TYPE  STRUCTURES 

A.   NUMERICAL  MODELING  OF  SURFACES  FOR  NEAR-FIELD 
PREDICTION 

There  are  two  principal  ways  in  which  the  method  of 
moments  can  be  used  to  model  antennas  together  with  the 
surrounding  structures.  The  simplest  way  to  model  a  solid 
surface  is  with  a  grid  of  wires,  the  so-called  wire-grid 
model.  The  other  approach  is  to  divide  the  surface  xnto 
patches  or  cells,  each  having  a  continuous  metallic  surface. 
This  is  the  so-called  surface-patch  model. 

The  successful  substitution  of  a  wire  grid  or  surface 
patch  for  a  continuous  metallic  surface  depends  upon  the 
fact  that  as  the  model  (the  wire  grid  or  surface  patch) 
becomes  smaller  relative  to  a  wavelength,  it  supports  a 
current  distribution  which  approximates  that  of  the 
continuous  surface.  The  current  is  only  an  approximation  to 
the  actual  current,  however,  and  as  such  it  can  be  expected 
to  accurately  predict  the  far  fields  but  not  necessarily  the 
near  fields.  This  is  due  to  the  fact  that  the  wire-grid 
model  supports  an  evanescent  reactive  field  on  both  sides  of 
its  surface.  An  actual  continuous  conducting  surface  is  not 
capable  of  supporting  such  a  field.  For  surface  patches  the 
current  is  assumed  to  be  continuous  over  a  patch  which  is 
only  a  simple  approximation  to  the  current  and  allows  no 
smooth  variation  over  the  body  being  modeled. 

8 


The  version  of  the  Numerical  Electromagnetic  Code  that 
will  be  used  in  this  work  is  NEC-3,  although  it  has  been 
reported  by  Burke  [Refs.  3  and  4]  that  there  are  some 
inaccuracies  in  NEC.  In  the  first  model,  the  wire  grid  with 
interior  dipole,  there  are  some  frequencies  in  which  the 
wavelengths,  are  more  than  ten-times  larger  than  the  side 
dimension  of  the  cell  segments,  and  this  may  affect  results. 
To  increase  confidence  in  the  accuracy  of  the  results,  two 
antenna  parameters  are  calculated:  the  average  power  gain 
and  the  input  impedance.  The  average  power  gain  provides  a 
good  check  of  the  solution  accuracy  if  it  is  close  to  one 
for  free  space  or  close  to  two  for  a  perfectly  conducting 
ground  plane.  For  our  purposes  we  will  consider  average 
power  gain  within  ±10%  of  theoretical  to  be  adequate  for 
engineering  purpose.  Correct  average  power  gain  does  not 
insure  accurate  near-field  prediction.  Therefore,  input 
impedance,  which  is  in  general  complex,  is  also  checked. 
The  real  part  of  the  impedance  relates  to  radiated  far-field 
power  while  the  imaginary  portion  relates  to  the  balance  of 
reactive  energy  in  the  near  field. 

In  the  first  model  analyzed,  the  wire  grid  with  a 
dipole  inside,  there  are  no  measured  values  for  comparison. 
In  the  second  model,  the  wire-grid  box  with  the  monopole, 
the  measured  admittance  is  available.  The  next  section 
describes  these  two  models. 


B.   WIRE  GRID  BOX  WITH  INTERIOR  DIPOLE 

A  WIRE-GRID  box  with  an  interior  dipole  was  modeled  to 
obtain  numerical  values  of  average  power  gain  and  input 
impedance.  The  model  constructed  was  a  2  meter  per  side 
cube-shaped  box,  with  a  dipole  inside,  1  Meter  long.  The 
grid  was  .4  by  .4  meter  cells,  with  each  cell  side  divided 
into  two  segments  (see  Figure  3.1).  (Initial  wire  grid  box 
models  employed  both  one  and  two  segments  per  cell  side  and 
favored  2  segments  per  side  as  opposed  to  the  conventional  1 
segment  per  side.)  The  model  was  run  at  selected  frequen- 
cies from  3  to  4  50  MHz,  in  free  space.  These  parameters  were 
chosen  to  span  geometry  parameters  where  guidelines  for  NEC 
wire  modeling  have  been  presented.   These  guidelines  were: 

1.  Ratio  between  segment  length  (D)  and  segment  radius 
(a)  greater  than  or  equal  to  two. 

2.  Angle  between  wires  greater  than  27r(D/a)~-'- 

3.  Spacing  between  segments  greater  than  or  equal  to  10a. 
Theoretically,  this  wire-grid  box  with  the  dipole  inside 

should  radiate.  Therefore,  average  power  gain  and  input 
impedance  will  be  calculated.  If  acceptable  values  are 
obtained,  a  statement  concerning  the  accuracy  of  near-field 
calculations  via  NEC  can  be  prepared.  In  this  case,  average 
power  gain  must  be  close  to  one  (free  space)  .  Values  of 
acceptable  input  impedance  must  be  limited  to  those  where 
the  ratio  of  the  real  to  the  imaginary  part  are  less  than 
10"*.  This  limit  is  because  of  the  accuracy  of  calculating 
complex  numbers  in  NEC. 
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1  MT.  DIPOLE  IN  2MT.  WIRE  GRID  BOX 


THETA  =  60.00    PHI  =  50.00    ETA  =  90.00 


Figure    3.1      1   Meter   Dipole   in   2   Meter  Wire- 
Grid   Box. 
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C.   MONOPOLE  ON  WIRE  GRID  BOX 

Bhattacharya,  Long  and  Wilton  [Ref.  2]  performed  an 
experimental  investigation  to  determine  the  input  admittance 
characteristics  of  a  monopole  antenna  mounted  on  a  con- 
ducting cubical  box  over  a  ground  plane.  In  this  case,  a  6 
cm.  monopole  was  positioned  at  the  center,  edge,  and  corner 
on  the  top  of  the  box  and  the  admittance  was  measured  (see 
Table  1).  This  measured  data  [Ref.  2]  was  compared  with 
calculated  data  obtained  using  a  five-sided  WIRE-GRID  BOX 
model  of  .1  meter  per  side  and  cells  of  .0125  by  .0125 
meters. 


Monopole 


Figure  3.2   Geometry  of  a  Monopole  Mounted  on  a  Cube  Shaped 
Box  Over  a  Ground  Plane.  [From  Ref.  2] 

12 


TABLE  1 
MEASURED  DATA  FOR  A  6  CM.  MONOPOLE  ON  A  CUBE  BOX  [REF.2] 


FREQUENCY 
GHZ. 

6  cm.  MONOPOLE  POSITION 

AT  CENTER 

AT  EDGE  S.SOcm 
FROM  CENTER 

AT  CORl^ER  5. 15 cm. 
FORM  CENTER 

ADMITTANCE 

1 .  0   GHZ  . 

5.00+J20.0 

7. 00+ j 17.0 

7.00+J14.0 

1.05  GHZ. 

14.0+J28.0 

1.05+J20.0 

12.0+J13.0 

1 . 1   GHZ  . 

37.0+J18.U 

24.0+J17.0 

18.0+jlO.O 

1.15  GHZ. 

41.0-j8.00 

32.0+J2.00 

21.0+J4.00 

1 .  2   GHZ  . 

26.0-J16.0 

26.0-j7.00 

20.0-j2.00 

1.225GHZ. 

19.0-J15.0 

26.0-jll.O 

17.0-j4.50 

1 . 3   GHZ . 

11.0-jlO.O 

12.0-j9.00 

12.0-j5.00 

1.4  GHZ. 

7.00-J6.00 

7.00-J5.00 

8.00-J3.00 

In  order  to  simulate  the  cube  box  in  Figure  3.2,  a  model 
of  a  6  cm.  monopole  divided  into  5  segments  was  placed  on 
top  of  the  wire-grid  box  at  the  center  (see  Figure  3.3),  at 
the  edge  3.75  cm.  from  the  center  (see  Figure  3.4),  and  at 
the  corner  5.3  cm.  from  the  center  (see  Figure  3.5).  The 
feed  of  the  antenna  was  placed  at  the  base  and  second 
segment  up  on  the  monopole  at  the  center  of  the  wire  grid 
only  to  see  if  there  was  any  major  difference  in  the 
results. 

The  calculated  average  power  gain  and  input  impedance 
are  shown  in  Tables  2 ,  3 ,  and  4 . 
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D.   MONOPOLE  ON  SURFACE  PATCH  BOX 

In  the  same  fashion  as  the  wire  grid  box,  a  five  sided 
SURFACE-PATCH  BOX  of  .1  meter  per  side  and  cells  .02  by  .02 
meters  was  constructed  in  order  to  simirlate  the  cubical  box 
described  in  Reference  1.  For  this  model,  two  different 
tops  of  the  cube  were  designed  to  allow  placement  of  the 
model  of  a  6  cm.  monopole  divided  into  five  segments  at  the 
center  (see  Figure  3.6),  at  the  edge  3.5  cm.  from  the  center 
(see  Figure  3.7)  and  at  the  corner  4.93  cm.  from  center  (see 
Figure  3.8).  In  each  case,  the  monopole  was  fed  at  the  base 
segment. 

The  calculated  average  power  gain  and  input  impedance 
are  shown  in  Tables  5,  6,  and  7,  and  description  of  the 
results  are  discuss  in  Chapter  V. 
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6  CM.  MONOPOLE  AT  CENTER  OF  WIRE  GRID  BOX 


THETA  =  60.00    PHi  =  60.00    ETA  =  SO. 00 


Figure    3.3       6    cm.    Monopole    at   Center    of   Wire- 
Grid   Box. 
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6  CM.  MONOPOLE  AT  EDGE  3.75  CM.  FROM  CENTER 


THETA  =  60.00    PHI  =  60.00    ETA  =  SO.OO 


Figure    3.4       6    cm.    Monopole    at    Edge    3.75    cm, 
from   Center. 
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6  Ch\.  MONOPOLE  AT  CORNER  5.3  C^.  FROM  CENTER 


THETA  =  70.00     FHI  =  55.00    ETA  =  90.00 


Figure  3.5   6  cm.  Monopole  at  Corner  5.3  cm, 
from  Center. 
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TABLE  2 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  CENTER 
OF  A  WIRE-GRID  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 

AVG . Power 
Gain 

IMPEDANCE   / 

ADMITTANCE 

1  Seg. 

2  Seg. 

1  Seg. 

2  Seg. 

AVG.l-2Seg. 

1.0 

GHZ. 

1.76 

1.87 

13.0-J40.7 
7.15+J22.2 

15.1-J45.4 
6.49+J19.6 

14.1-J43.0 
6.82+J20.9 

1.05 

GHZ. 

1.79 

1.91 

16.8-j23o3 
20.3+J28.2 

17.9-J25.5 
18.4+J26.3 

17.4-J24.4 
19.4+J27.3 

1.1 

GHZ. 

1.8 

1.93 

21.6-j6.60 
42.4+J13.0 

24.l-j7.00 
38.2+jll.l 

22.9-j6.80 
40.3+J12.0 

1.15 

GHZ. 

1.82 

1.94 

27.4+J9.30 
32.7-jll.l 

29.2+J10.3 
30.4-J10.7 

28.3+J9.80 
31.6-J10.9 

1.2 

GHZ. 

1.83 

1.96 

34.3+J24.1 
19. 5- j 13. 7 

36.9+J27.4 
17.4-J13.0 

35.6+J25.8 
18.5-J13.4 

1.3 

GHZ. 

1.84 

1.97 

50.3+J50.0 
lO.O-jlO.8 

52.1+J56.2 
8.88-J9.57 

51.2+J53.1 
9.44-J10.2 

1.4 

GHZ. 

1.86 

1.98 

66.7+J70.6 
7.07-J7.48 

66.8+J79.0 
6.24-J7.38 

66.8-f-j74.8 
6.70-J7.43 
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TABLE  3 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  EDGE 
OF  A  WIRE-GRID  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 

AV. POWER 
GAIN 

IMPEDANCE/ 
ADMITTANCE 

1 . 0    GHZ . 

1.-17 

28.9-J45.8 
9.90+J15.6 

1.05   GHZ. 

1.10 

31.4-J30.8 
16.2+J16.0 

1 . 1    GHZ . 

1.08 

34.4-J16.0 
24.0+jll.l 

1.15   GHZ. 

1.11 

38.l-jl.08 
26.2+J0.74 

1.2    GHZ. 

1.18 

43.0+J13.5 
21.2-j6.67 

1.3    GHZ. 

1.40 

55.4+J41.0 
ll.7-j8.63 

1 .  4    GHZ . 

1.65 

71.1+J65.7 
7.60-J7.00 
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TABLE  4 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  CORNER 
OF  A  WIRE-GRID  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 

AVG.  POWER 
GAIN 

IMPEDANCE/ 
ADMITTANCE 

1 . 0    GHZ . 

1.86 

48.0-J47.4 
10.6+J10.4 

1.05   GHZ. 

1.86 

49.9-J34.5 
13.7+J9.37 

X . 1    GHZ . 

1.87 

52.0-j6.63 
16.6+J6.61 

1.15   GHZ. 

1.87 

55.3-j6.63 
17.8+J2.13 

1 . 2    GHZ . 

1.87 

60.1+J7.60 
16.4-j2.07 

1 . 3    GHZ . 

1.88 

74.3+J34.6 
ll.l-j5.10 

1 . 4    GHZ . 

1.88 

92.7+J55.0 
8.00-J4.73 
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6  CM.  MONOPOLE  AT  CENTER  OF  SURFACE  PATCH  BOX 


THETA  =  60.00  FHi  =  60.00  ETA  =  90.00 


Figure  3.6   6  cm.  Monopole  at  Center  of  Surface- 
Patch  Box. 
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6  CM.  MONOPOLE  AT  EDGE  3.5  CM.  FROM  CENTER 


THETA  =  60.00    PHI  =  60.00    ETA  =  90.00 


Figure    3.7      6   cm.    Monopole   at   Edge   3.5   cm, 
from   Center. 
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6  CM.  MONOPOLE  AT  C0RNER4.93  CM.  FROM  CENTER 


THETA  =  70.00     PHI  =  55.00    ETA  =  SO. 00 


Figure   3.8      6   cm.    Monopole   at   Corner   4.93    cm, 
from   Center. 
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TABLE  5 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  CENTER 
OF  A  SURFACE-PATCH  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 


AVG. 
GAIN 


POWER 


IMPEDANCE/ 
ADMITTANCE 


1.00   GHZ 


2.02 


10.8- 
20. 6+ j 


20.2 
38.4 


1.05   GHZ 


2.01 


1.10   GHZ 


2.01 


14.3-_ 
60. 6+ j 


5.49 
23.2 


18. 8+ j 
43.9- 


8.58 
19.9 


1.15   GHZ 


2.01 


24.2  + 
22.7- 


21.9 
20.0 


1.175  GHZ 


2.01 


27.3  + 
17.1- 


28.2 
18.2 


1 . 2    GHZ 


2.01 


30.6+j 
14.5- 


34.2 
16.1 


1.225  GHZ 


2.01 


34.0+ 
12.3- 


40.0 
14.5 


1.30   GHZ 


2.01 


44.8+ 
8.84- 


55.5 
10.9 


1.4    GHZ 


2.02 


58.4+J71.8 
6.81-J8.37 
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TABLE  6 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  EDGE 
OF  A  SURFACE-PATCH  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 

AVG.  POWER 
GAIN 

IMPEDANCE/ 
ADMITTANCE 

1.00   GHZ. 

1.92 

22.6+J7.40 
39.0-J13.1 

1.05   GHZ. 

1.92 

25.3+J18.4 
25.9-J18.8 

1.10   GHZ. 

1.92 

28.3+J29.4 
16.9-J17.5 

1.15   GHZ. 

1.88 

32.9+J40.0 
12.1-J14.8 

1.175  GHZ. 

1.88 

35.2+J46.0 
10.4-J13.7 

1 . 2    GHZ . 

1.88 

37.6+J51.5 
9.23-J12.6 

1.225  GHZ. 

1.88 

40.3+J57.0 
8.77-J11.7 

1 . 3    GHZ . 

1.89 

49.5+J72.3 
6.44-J9.41 

1.4    GHZ. 

1.89 

123.+J220. 
1.94-J3.46 
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TABLE  7 


CALCULATED  DATA  FOR  A  MONOPOLE  AT  THE  CORNER 
OF  A  SURFACE  PATCH  BOX,  BASE  SEGMENT  FEED 


FREQUENCY 
GHZ. 

AV. POWER 
GAIN 

IMPEDANCE/ 
ADMITTANCE 

1.00   GHZ. 

1.85 

34.1-J53.2 
8.52+J13.3 

1.05   GHZ. 

1.86 

36.2-J40.5 
12.2+J13.7 

1.10   GHZ. 

1.87 

38.7-J20.4 
17.0+J12.1 

1.15   GHZ. 

1.83 

43.0-J14.2 
20.9+J6.93 

1.175  GHZ. 

1.83 

45.l-j7.38 
21.5+J3.52 

1.2  0   GHZ. 

1.83 

47.6-j .622 
21.0+J.270 

1.225  GHZ. 

1.83 

50.4+J6.01 
19.5-j2.33 

1 . 3    GHZ . 

1.83 

60.6+J24.7 
14.l-j5.70 

1 . 4    GHZ . 

1.82 

77.4+J45.9 
9.54-J5.66 
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IV.   WIRE  GRID  SOURCE  BOX 

The  results  of  the  average  power  gain  and  input 
impedance  of  the  wire-grid  cube-shaped  box  model  are 
displayed  in  Table  8.  For  all  of  the  frequencies  chosen, 
the  average  power  gain  was  close  to  one.  (This  model  was  in 
free  space  where  the  theoretical  average  power  gain  must  be 
one)  .  Only  input  impedance  values  for  frequencies  from  42 
to  450  MHz  can  be  considered  useable  because  the  ratios  of 
the  real  part  to  the  imaginary  part  are  less  than  10"^.  For 
frequencies  of  150  and  300  MHz,  where  the  side  dimension  of 
the  box  is  one  and  one-half  wavelength,  respectively,  the 
ratio  of  the  real  part  to  the  imaginary  part  was  10^,  which 
indicates  that  the  dipole  is  resonant  (A/2  and  /\  long)  and 
radiates  very  well  through  the  widely-spaced  wire  grids. 
However  when  the  model  gets  smaller  in  relation  to  the 
wavelength,  as  for  30  MHz  and  below,  the  results  have  ratios 
of  lO'^  and  higher  between  the  real  and  the  imaginary  parts 
of  the  input  impedance. 

Calculations  of  the  electrical  near  field  for  the 
frequencies  of  30,  75  and  300  MHz,  were  performed  for  a  path 
starting  at  the  interior  dipole  and  passing  through  the  wire 
grid  to  the  outside  of  the  box,  (see  Figures  4.1  and  4.2  and 
Appendix  A)  .  The  plots  show  the  E  field  along  the  Y  and  Z 
axis.    It  can  be  seen  in  all  the  plots  that  the  fields 
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inside  the  box  at  the  wire  grid  walls  are  down  several 
orders  of  magnitude  from  those  at  the  dipole's  surface. 
Passing  through  the  grid,  the  fields  drop  suddenly  by 
another  order  of  magnitude  for  the  two  lower  frequencies.  At 
300  MHz,  where  the  grid  spacing  is  very  wide  (.4A)  the  box 
is  quite  transparent,  with  little  shielding  effect. 
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TABLE  8 


CALCULATED  DATA  FOR  THE  DIPOLE  INSIDE 
THE  WIRE-GRID  BOX 


FREQUENCY 
Mhz. 

LAMBDA 
Mts. 

AV.  POWER 
GAIN 

IMPEDANCE 

3 

100 

.98 

.88 

10^  -j.22 

105 

7.5 

40 

.98 

.53 

10_2-j .90 

104 

15 

20 

.98 

.18 

lO'l-j .45 

104 

30 

10 

.99 

.41 

10"2-j .21 

104 

35 

8.57 

.97 

.11 

10"l-j .18 

104 

42 

7.14 

.97 

.22 

10°  -j.l5 

104 

50 

6 

.98 

.99 

10°  -j.l2 

104 

62 

4.83 

.98 

.47 

10^  -j .89 

103 

75 

4 

.98 

.33 

102  -j.62 

103 

92 

3.26 

.98 

.74 

102  -j.54 

103 

110 

2.72 

.98 

.19 

102  -j.30 

103 

150 

2 

.98 

.53 

102  +J.33 

102 

300 

1 

.97 

.18 

10^  +j .12 

104 

450 

.66 

.97 

.12 

10^  +J.60 

102 
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Figure    4.1      E-Field   along   Dipole   Axis    (Z-Axis)     for    IM. 
Dipole    inside   2M.      Wire-Grid   Box   at   30   MHz. 
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Figure  4.2  E-Field  perpendicular  to  Dipole  Axis  (Y-Axis) 
for  IM.   Dipole  inside  2M.   Wire  Grid  Box  at 
30  MHz. 
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V.   RESULTS  OF  THE  VARIATION  OF  THE  POSITION 
OF  THE  MONOPOLE  ON  THE  BOX 

A.   WIRE  GRID  CASE 

To  describe  the  curves  which  show  the  measurements  and 
the  calculation  of  the  input  impedance  for  the  variation  of 
the  position  of  the  6  cm.  monopole,  five  composite  graphs 
were  constructed  showing  the  curves  of  NEC-calculated  and 
measured  Conductance  (G)  and  Susceptance  (B)  from  Reference 
2,  (see  Figures  5.1  through  5.5).  In  the  case  of  the 
monopole  at  the  center  of  the  box,  the  resonant  frequency, 
as  defined  by  the  axis  crossing  of  the  susceptance  curve, 
decreased  from  1.138  to  1.125  GHz  for  Figure  5.1  (base 
segment  feed),  decreased  from  1.138  to  1.122  GHz  for  Figure 
5.2  (second  segment  feed),  and  decreased  from  1.138  to  1.123 
GHz  for  Figure  5.3  (the  average  of  the  two  different  feed 
methods) .  All  these  are  relative  to  the  measured  suscep- 
tance from  Reference  2.  The  variation  was  small  due  to  the 
feed  position.  The  feed  segment  was  varied  only  for  this 
position  of  the  monopole.  In  some  cases,  when  the  feeding 
segment  is  tied  to  many  other  segments,  errors  may  result. 
When  feeding  a  vertical  segment  tied  to  four  horizontal 
segments,  calculations  are  difficult  because  the  segments 
which  are  fed  usually  like  to  have  a  straight  segment  on 
both  sides.  Due  to  the  good  results  obtained  with  the  base 
segment   feed,   all   the   other  comparisons   are  with   the 
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monopole  base  segment  feed;  the  measured  data  was  also 
obtained  by  feeding  that  segment.  The  peak  of  the  real  part 
of  the  input  admittance  was  lower  for  the  calculated  value, 
changing  from  4  6.4  to  41  mS. 

In  the  case  of  the  monopole  at  the  edge  of  the  box,  the 
calculated  resonant  frequency  decreased  slightly  from  1.164 
to  1.155  GHz  and  the  peak  conductance  varied  from  31.7  to 
26.2  mS  (Figure  5.4)  in  comparison  to  measurements.  In  the 
case  of  the  monopole  at  the  corner,  the  NEC  resonant 
frequency  decreased  from  1.192  to  1.157  GHz  and  the  peak 
conductance  from  2  0.9  to  17.8  mS  (Figure  5.5).  For  this 
wire-grid  model,  the  results  of  the  calculated  data  are  very 
close  to  the  measured  values  from  Reference  2.  Also,  the 
result  of  the  average  power  gain  for  this  model  were  close 
to  two  (correct  value  of  average  power  gain  for  perfect 
conducting  ground  plane) . 

B.   SURFACE-PATCH  CASE 

For  this  model,  the  results  of  the  calculated  (NEC) 
input  admittance  were  very  good  for  the  monopole  located  at 
the  corner.  There  was  a  shift  in  resonant  frequency  of 
about  -5  %  for  the  monopole  at  center,  (see  Figure  5.6).  For 
the  edge-mounted  geometry,  the  correlation  between 
measurement  and  calculation  was  quite  poor  (see  figure  5.7). 
The  calculated  values  of  the  average  power  gain  were  close 
to  two  for  all  positions  of  the  monopole  (see  Tables  5,  6, 
and  7,  Chapter  III).   For  the  monopole  at  the  center  of  the 
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surface-patch  box,  the  resonant  frequency  decreased  from 
1.138  to  1.052  GHz  and  the  peak  conductance  increased  from 
46.4  to  60.6  mS  (see  Figure  5.6).  The  resonant  frequency 
for  corner  mounting  changed  from  1.192  to  1.2  GHz  and  the 
peak  conductance  increased  from  20.9  to  21.5  mS  (see  Figure 
5.8).  This  accuracy  might  be  due  to  the  fact  that  in  this 
case  the  top  of  the  surface-patch  box  had  been  more  finely 
divided  into  smaller  patches  than  for  the  other  two  mounting 
cases.  This  was  done  in  order  to  position  the  monopole  at  a 
specific  distance  from  the  center,  close  to  the  distance 
used  in  Reference  2. 
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Figure  5.1 


6  cm.  Monopole  at  center  of  Wire-Grid 
Box,  Base  Segment  Feed. 
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Figure  5.2   6  cm.  Monopole  at  center  of  Wire-Grid 
Box,  Second  Segment  Feed. 
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Figure  5.3   6  cm.  Monopole  at  center  of  Wire-Grid  Box, 
Average  of  Base  and  second  Segment  Feeds. 
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Figure  5.4 


6  cm.  Monopole  at  Edge  of  Wire-Grid 
Box,  Base  Segment  Feed. 
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Figure  5.5 


6  cm.  Monopole  at  Corner  of  Wire-Grid 
Box,  Base  Segment  Feed. 


39 


X 

o 

n: 

LL 

o 

01 
H 
Z 
CJ 


o 

CL 

o 
o 


LLi 


o 


09 


Ot'  OC  OS  01  0  01- 

(SOHl/Nmil/^)  33NVllllAiaV 


Figure    5.6      6    cm.    Monopole    at   Center   of   Surface-Patch 
Box,    Base   Segment   Feed. 
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Figure  5.7 


6  cm.  Monopole  at  Edge  of  Surface-Patch 
Box,  Base  Segment  Feed. 
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Figure  5.8 


6  cm.  Monopole  at  Corner  of  Surface-Patch 
Box,  Base  Segment  Feed. 
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VI.   CONCLUSION  AND  RECOMMENDATIONS 

The  initial  goal  of  this  work  was  to  investigate  the 
accuracy  of  near-field  calculations  using  the  Numerical 
Electromagnetic  Code  (NEC)  on  ship-like  structures.  Since  no 
validation  bench  mark  results  were  available  for  ship 
antenna  installations,  a  modeling  exercise  was  undertaken, 
which  would  yield  insight  for  near-field  work  in  the  future. 
Box-like  structures  were  analyzed,  one  with  an  interior 
dipole  and  one  with  a  monopole  mounted  on  the  top.  For  the 
dipole  inside  the  wire  grid,  there  were  no  measured  values 
to  compare  to,  but  the  calculations  of  the  average  power 
gain  and  input  impedance  were  as  expected.  The  E  field 
variation  was  also  examined  as  a  function  of  frequency.  The 
fields  inside  the  box  drop  off  as  the  observation  point 
moves  away  from  the  source.  For  a  non-radiating  structure, 
the  fields  would  be  similar  to  the  interior  of  a  small 
cavity,  with  standing  wave  patterns.  Any  input  resistance 
in  the  perfect  cavity  case  could  only  come  from  dipole 
conductor  loss.  The  leaky  nature  of  the  wire  box  produced 
non-zero  input  resistance  for  the  dipole,  but  as  the 
frequency  was  lowered  to  a  point  where  the  grid  size 
approached  ^/20,  the  real  part  of  the  input  impedance  was 
numerically  swamped  out  by  the  reactance,  an  expected  result 
arising  from  less  radiation  leakage.   The  value  of  the  leaky 
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box  model  is  that  it  raises  concerns  about  wire  grids  used 
as  models  of  ships  structures  for  HF  RZUDHAZ  predictions.  On 
a  real-world  ship,  the  antennas/sources  will  be  located  on 
the  exterior  of  the  ship.  The  predicted  near  fields  around 
these  radiators  may  be  subjected  to  substantial  calculation 
errors  if  the  wire  grid  surfaces  are  not  "tight"  enough. 
Measurements  of  near-fields  for  well-controlled  generic 
shipboard  shapes  are  needed  as  a  validation  benchmark  for 
additional  numerical  model  developments. 

In  the  case  of  the  monopole  over  the  cube-shaped  box, 
the  two  modeling  methods  that  NEC  has  available  were  tried. 
The  wire-grid  model  worked  very  well,  but  the  surface  patch 
gave  good  results  only  for  the  monopole  positioned  at  the 
corner  of  the  box  where  the  patches  at  the  top  are  smaller. 
If  smaller  patches  (finer  surface  current  definition)  for 
the  other  two  positions  of  the  monopole  had  been  used  the 
results  might  have  been  better.  Those  data  sets  are  easy  to 
produce,  but  will  be  difficult  to  calculate  because  of 
excessive  computer  time  requirements.  It  is  recommended 
that  this  be  done  on  a  super-computer  in  the  continuation  of 
this  work. 

This  study  is  an  important  step  in  the  direction  of 
near-field  modeling  of  shipboard  shapes  because  it  shows  the 
kind  of  wire-grid  definition  needed  to  get  accurate  input 
impedance  for  antennas  on  boxes  and  gives  an  idea  of  how  to 
apply  surface-patch  modeling  for  accurate  input  impedance. 
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APPENDIX   A 

PLOTS    OF    E-FIELDS    ALONG    DIPOLE    AXIS     (Z-AXIS)     AND 
PERPENDICULAR   AXIS     (Y-AXIS)     FOR    IM    DIPOLE    INSIDE 

WIRE    GRID    BOX 
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Figure  A.l   E-Field  Along  Dipole  Axis  (Z-Axis)  for  IM. 
Dipole  Inside  2M.  Wire  Grid  Box  at  75  MHz. 
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Figure  A. 2   E-Field  Perpendicular  to  Dipole  Axis  (Y-Axis) 
For  IM.  Dipole  Inside  2M.  Wire  Grid  Box  at  75  MHz. 
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Figure   A. 3 


E-Field  Along  Dipole  Axis  (Z-Axis)  for  IM. 
Dipole  Inside  2M.  Wire  Grid  Box  at  300  MHz 
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Figure  A. 4   E-Field  Perpendicular  to  Dipole  Axis  (Y-Axis) 
for  IM.  Dipole  Inside  2M.  Wire  Grid  Box  at  300  MHz. 
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APPENDIX  B 
INPUT  DATA  FILES 


CM  WIRE  GRID  10  CM.  LONG 

CM  CELLS  . 125  X. 125  GROUND 

CE 

GW  100,4,. 5,-. 5,1,0,-. 5,1..  01 

GM  1,7,0,0,0,0,0,-. 125,100. 100 

GW  260,8,. ^,-. ^,i.. 5,-. 5,0,. 01 

GM  2,3,0,0,0,0, . 125.0,200. 201 

GM  3,3,0,0,0,-. 125,0,0,200. 201 

GW  300,4, . 5,-. 5.1, . 5.0, 1,. 01 

GM  4,7,0,0,6,0.6,-. 125,300. 303 

GW  460,4,. 5,-. 375,1,0,-. 3 75,1,. 01 

GM  5,2,0,0,6,0,. Iz5, 0,400. 464 

GW  560,4,. 575, -. 5,1,. 375,0, 1,. 01 

GM  6,2,0,0,0,-. 120,6,0,560. 565 

GX  0, 110 

GS  0,0, . 1 

GE    1 

FR  0,0,0,0,1000 

GN  1 

WG 

XO  0 

EN 


1.   Data  File  for  Wire-Grid  Box. 


CM  6  Cr:.  MONOPOLS  AT  CENTER  OF  WIRE  GRID  FREO  1.  0  GHZ 
CM  CALCUL.'.TE  OF  AVERAGE  GAIIJ  AMD  IMPUT  IMPEDANCE 


o  V 


GW  2,8,G,-. 5,1,0^. -,!,. 01 

GW  3,8,. ^01-. 0,0,1,. 01 

GW  1, 5,0,6, 1,6,0, 1.  6, . 016 

G3  0,0, . 1 

EX  0,1,1,0,1,0 

PP  0,31, 4, 1062,0,0,3, 15 

::o  0 

EK 


2.   Data  File  6  cm.  Monopole  at  Center  Wire  Grid  Box. 
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CM  6  CM.  MONOPOLE  AT  EDGE  OF  WIRE  GRID  FREO  1.  0  GHZ 

CM  CALCULATE  OF  AVERAGE  GAIN  AND  IMPUT  IMPEDANCE 

CE 

GF 

GW  2,8,0,-.  5,1,0, . 5, 1, .  01 

GW  3,8,  .  5.0,  1,-.  5,0.1,.  01 

GW  1,5,-.  3 7d, 6,1,-. 375,0, 1. 6, . 016 

G3  0,0,  .  1 

EX. 0,1. 1,0,  1,0 

RP  0,31, 13, 1602,0,0,3, 15 

XO  0 

Ef] 


3.   Data  File  6  cm.  Monopole  at  Edge  Wire-Grid  Box. 


CM  5  CM.  MONOPOLE  AT  CORNER  OF  WIRE  GRID  FREO  1.0  GH: 

CM  CALCULATE  Or  AVERAGE  GAIN  AND  IMPUT  IMPED/^NCE 

CE 

GF 

GW  2,3,0^-. 5, 1,0, .5,1, .01 

GW  3,8,. D,0,i.-.^.Q,1,. 01 

C-vl    1,1;,.  57 Z-,  .  j75,1,.  375,.  375,1.  6,.  016 

G3  0,0, .  1 

GE 

EX  0  1  "I  0  "■  0 

:■-?  0'3i7i3ii602,0,45,3,  15 

XO  0 

-N 


4.  Data  File  6  cm.  Monopole  at  Corner-Wire  Grid  Box, 
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CE  SURFACE  PATCH   GROUND   FREQ.  1.  OGHZ 

SM3,3,  .5, -.5,0,  .5,-. 5,1 

SC  0,0,  0,-. 5,1 

SM  3,3,  . D,-. 5,0,  .5,0,0 

SC  0,0,  . 5,0,1 

OX   0,1 io 

SM  5,5,  -.5, -.5,1,  .5,-. 5,1 
SC  0,0,  . 5,. 5,1 
GS  0,0,. 1 

GE  1 

FR  0,0,0,0,1000 

GN  1 

v;g 

XO  0 


5.  Data  File  Surface-Patch  Box  for 
Monopole  at  Center. 


CE  SURFACE  PATCH   GROUND 

SM  3,3,  .5,-. 5,0,  .5, -.5,1 

SC  0,0,  0,-. 5,1 

SM  3,3,  .5,-.^,0,  .5,0,0 

SC  0,0,  . 5,0,1 

GZ  0, lio 

SM  5,5,  -.5, .5,1,  -.5, -.5,: 

SC  0,0,  0,-. 5, 1 

SM  10,^,  6, . 5,1,  0,-. 5, 1 

SC  0,0,  .  5, -.  5, 1 

GS  0 , 0 , . 1 

GE  1 

FR  0,0,0,0, 1000 

GN  1 

V/G 

XO  0 

EN 


Data  File  Surface-Patch  Box  for 
Monopole  at  Edge  and  Corner. 
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CM  6  CM.  MONOPOLE  AT  CENTER  SURFACE  PATCH   FREO  1.  0  GHZ 

CM  CALCULATE  Or  AVERAGE  GAIN  AND  IMPUT  IMPEDACE. 

CE 

GF 

GW  1,5,0,0,1,0,0,1. 6,. 016 

GS  0 , 0 , . 1 

GS 

EX  0, 1,1,0,1.0 

R?  0,31,4,1002,0,0,3,15 

XO  0 

E5I 


7.   Data  File  for  6  cm.  Monopole  at 
Center  Surface-Patch  Box. 


CM  6  CM.  MONOPOLE  AT  EDGE  SURFACE  PATCH   FREQ  1.  0  GHZ 

CM  CALCULATE  OF  AVERAGE  GAIN  AND  IMPUT  IMP ED.nCE. 

CE 

GF 

GW  1,5,-, 35,0,1,-. 35,0,1. 6,. 016 

GS  0,0, . 1 

GE 

EX  0  "I  ■>  0  1  0 

R?  0!3i, 13^1602,0,0,3, 15 

XO  0 

—0 


8.   Data  File  for  6  cm.  Monopole  at 
Edge  Surface  Patch. 


CM  MONOPOLE  6CM.  AT  CORNER  SURFACE  PATCH   FREO  1.  0  GHZ 

CM  CALCULATE  OF  AVERAGE  GAIN  AND  INPUT  IMPEDACE. 

CE 

GW  1,5, . 3  5,. 3  5, 1, . 35, . 5  5, 1. 6, . 016 

GS  0, 0, .  1 

GE 

EX  0, 1,1,0, 1,0 

R?  0,3i, i3, 1602,0,45,3, 15 

XO  0 

EN 


9.   Data  File  for  6  cm.  Monopole  at 
Corner  Surface  Patch. 
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CM  AVERAGE  GAIN  DIPOLE  IN  A  WIRE  GRID 

CM  CELLS  . 4X. 4  DIPOLE  LENGTH  IZ  FREE  SPACE 

CE 

GW  100,1,1,-1,1,0,-1,1..  04 

GM  1,2,0,0,0,0,. 4,0, 106.  100 

GH  1,2,0,0,0,0,0,-. 4,100.  100 

GV7  200,  1,1,-1,1,1.0.  1,.  04 

GM  2,2,0,0,0, -. 4,0,6,200.  200 

GM  2,2,0,0,0,0,0,-.  4,200. 200 

GW  300, 1, 1,-1,1, 1.-1,0,  .  04 

GM  3,2,0,0,0,-. 4,0,0,300.  300 

GM  3,2,0,0,0,0, . 4,0,300.  300 

GX  0 , 1 1 1 

GE 

FR  0,0,0,0,15 

WG 

EN 


CM  DIPOLE  IN  WIRE  GRID  AT  0,0,-.  5   0,0,.  5 

Ci^I  FREQ.  15  MHZ,  CALCULATE  AVERAGE  GAIN,  INPUT  IMPEDANCE 


GW  1,5,0,0,-. 5,0,0, . 5, . 001 

GE 

^/  0  ""  '^0  1  0 

RP  o;3i,4, ic62,0,0,3, 15 

XO  0 

EN 


10.   Data  Files  Wire-Grid  Box  2  by  2  Meters 
Dipole  1  Meters  Long. 
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APPENDIX  C 
DESCRIPTION  OF  NEC 
The  Numerical  Electromagnetics  Code  (NEC)  is  a  user- 
oriented  computer  code  for  analysis  of  the  electromagnetic 
response  of  antennas  and  other  metal  structures.  It  is 
built  around  the  numerical  solution  of  integral  equations 
for  the  currents  induced  on  the  structure  by  sources  or 
incident  fields.  This  approach  avoids  many  of  the  simpli- 
fying assumptions  required  by  other  solution  methods  and 
provides  a  highly  accurate  and  versatile  tool  for  electro- 
magnetic analysis. 

The  code  combines  an  integral  equation  for  smooth 
surfaces  with  one  specialized  to  wires  to  provide  for 
convenient  and  accurate  modeling  of  a  wide  range  of 
structures.  A  model  may  include  nonradiating  networks  and 
transmission  lines  connecting  parts  of  the  structure, 
perfect  or  imperfect  conductors,  and  lumped  element  loading. 
A  structure  may  also  be  modeled  over  a  ground  plane  that  may 
be  either  a  perfect  or  imperfect  conductor. 

The  excitation  may  be  either  voltage  sources  on  the 
structure  or  an  incident  plane  wave  of  linear  or  elliptic 
polarization.  The  output  may  include  induced  currents  and 
charges,  near  electric  or  magnetic  fields,  and  radiated 
fields.  Hence,  the  program  is  suited  to  either  antenna 
analysis  or  scattering  and  EMP  studies. 
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The  integral  equation  approach  is  best  suited  to  struc- 
tures with  dimensions  up  to  several  wavelengths.  Although 
there  is  no  theoretical  size  limit,  the  numerical  solution 
requires  a  matrix  equation  of  increasing  order  as  the 
structure  size  is  increased  relative  to  the  wavelength. 
Hence,  modeling  very  large  structures  may  require  more 
computer  time  and  file  storage  than  is  practical  on  a 
particular  machine.  In  such  cases  standard  high-frequency 
approximations  such  as  geometrical  optics,  physical  optics, 
or  geometrical  theory  of  defraction  may  be  more  suitable 
than  the  integral  equation  approach  used  in  NEC. 

NEC  contains  the  Numerical  Green's  Function  for 
partitioned-matrix  solution  and  a  treatment  for  lossy 
grounds  that  is  accurate  for  antennas  very  close  to  the 
ground  surface.  It  also  includes  an  option  to  compute 
maximum  coupling  between  antennas  and  useful  options  for 
structure  input . [ Ref . 6 : p . 1 ] 
1.   ZONING  CONSIDERATIONS 

NEC  is  a  discrete  sampling  code  where  a  complex 
structure  must  be  discretized  into  a  number  of  simple 
elements  (wires  or  surface  patches)  to  which  the  Electrical 
Field  Integral  Equation  (EFIE)  or  Magnetic  Field  Integral 
Equation  (MFIE)  are  applied.  As  with  any  description  of  the 
real  world,  there  are  approximations,  but  because  of  the 
versatility   in  modeling  the   geometry  of  a   structure, 
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its  approximations  more  closely  resemble  nature.  The 
resemblance  is  strongly  influenced  by  the  choice  of  zoning 
(i.e.,  dissecting)  the  structure  in  the  program.  The 
smaller  the  geometric  elements,  the  closer  the  model  comes 
to  reality.  However,  the  smaller  the  elements,  the  larger 
number  of  elements,  which  means  the  larger  the  matrix  of 
equations  and  hence,  the  more  costly  the  solution.  There  is 
a  point  beyond  which  smaller  zones  do  not  yield  a 
substantially  more  accurate  solution;  it  may  be  necessary  to 
rerun  the  problem  with  increasingly  smaller  elements  to  find 
the  point  of  diminishing  returns.  The  choice  of  proper 
zoning  then  is  gained  with  experience  and  becomes  as  much  of 
an  art  as  it  is  a  science.  The  guidelines  for  the  science 
aspect  are  as  follows. [Ref. 6] 
a.    Wires 

Segments  should  follow  the  paths  of  conductors 
using  a  piece-wise  linear  fit  on  curves.  Generally, 
segments  lengths  (^)  should  be  less  than  .  1  A  ;  shorter 
segments  (.05^  or  less  may  be  needed  at  critical  regions 
(junctions  or  curves)  .  Segments  smaller  than  10"-^  y\  should 
be  avoided  since  the  similarity  of  constant  and  cosine 
components  lead  to  numerical  inaccuracy.  The  radius  of  the 
wire  (a)  relative  to  ^  depends  on  the  Kernel  use  in  the 
Integral  Equation.  Two  options  exist.  The  thin-wire  Kernel 
models  a  filament  current,  while  the  extended  Kernel  models 
a  uniform  current  distribution  around  the  segment  surface. 
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The  field  of  the  distributed  current  is  approximated  by  the 
first  two  terms  in  a  series  expansion  of  the  exact  field,  in 
powers  of  a^ .  The  first  (a^)  term  is  identical  to  the  thin- 
wire  Kernel ;  the  second  term  extends  the  accuracy  for  larger 
values  of  a.  Both  Kernels  incorporate  the  thin-wire 
approximations  and  both  require  27ra/^«l.  The  thin-wire 
Kernel  requires  A  /a  >  8 ;  the  extended  Kernel  relaxes  this 
to  A /a  >  2.  These  values  ensure  errors  are  less  than  1% . 
The  extended  Kernel  is  used  at  free  wire  ends  and  between 
parallel,  connected  segments.  The  thin-wire  Kernel  is 
always  used  at  bends, 
b.     Surfaces 

A  conducting  surface  is  modeled  by  small  flat 
surface  patches  which  conform  as  closely  as  possible  to 
curved  surfaces.  The  parameter  defining  a  patch  is  a  normal 
unit  vector,  originating  from  the  center  of  the  patch, 
defined  in  Cartesian  coordinates.  Each  patch  must  be 
connected  by  a  wire  at  the  patch  center  for  the  program  to 
integrate  the  surface  current.  The  code  divides  each  patch 
into  four  equal  patches  about  the  wire  end,  along  the  unit 
vector  lines  describing  the  surface  of  the  of  the  patch.  An 
interpolation  function  is  applied  to  the  four  patches  to 
represent  placement  of  the  current  onto  the  patches,  and  the 
function  is  numerically  integrated.  Patches  with  wires 
connected  to  them  (active  patches)  should  be  chosen  to  be 
approximately  square  with  sides  parallel  to  the  unit  vectors 
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defining  the  surface.  Only  one  wire  may  connect  to  a  patch, 
that  wire  may  not  be  connected  to  another  patch.  A  minimum 
of  about  25  patches  should  be  used  per  wavelength  of  surface 
area;  the  maximum  size  for  an  individual  patch  is  about  .04 
square  wavelength.  The  number  of  patches  used  increases, 
and  the  size  of  each  patch  decreases,  as  the  radius  of 
curvature  decreases.  Smaller  patches  should  be  used  at 
edges  since  the  current  amplitude  may  vary  rapidly  in  this 
region.  Long  narrow  patches  should  be  avoided.  Patches  are 
restricted  to  modeling  voluminous  bodies  with  closed 
surfaces;  parallel  surfaces  on  opposite  sides  cannot  be  too 
close  together. 

c.     Ground  Plane 

For  a  perfectly  conducting  ground,  the  code  gen- 
erates a  reflected  image.  Structures  may  be  close  to,  or 
contact,  the  ground;  however,  for  a  horizontal 
wirezyh^  +a^  >10"^ A  where  a  =  wire  radius,  h  =  height  of  wire 
axis  above  the  ground  plane,  and  h/a  >  3. 

A  finitely  conducting  ground  may  be  modeled  by  an 
image  modified  by  the  Fresnel  Plane-wave  Reflection 
coefficients.  This  method  is  fast,  but  of  limited  accuracy 
and  should  not  be  used  for  structures  close  to  the  ground, 
or  having  a  large  horizontal  extent  over  the  ground.  The 
Sommerfeld/Norton  model  uses  the  exact  solution  and  is 
accurate  close  to  the  ground;  the  horizontal  restriction  is 
the  same  as  for  a  perfect  ground.   This  model  is  only  used 
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for  wire  to  wire  interactions,  for  surfaces  the  code  reverts 
to  Fresnel  Reflection  coefficients.  [Ref.  6]. 
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APPENDIX  D 
RADIATION  HAZARD  (RADHAZ)  ON  BOARD  SHIPS 
The   near   fields   of   Navy   antennas   are   primarily 
important  as  radiation  hazards  (RADHAZ) .   There  are  essen- 
tially three  types  of  radiation  hazards: 
RADHAZ  to  Personnel  (HERP) 
RADHAZ  to  Ordnance  (HERO) 
RADHAZ  to  Fuel  (HERE) 

These  hazards  and  their  general  near  field  criteria  are 
discussed  in  the  following  sections. 

a.   Hazards  of  Electromagnetic  Radiation  to 
Personnel  (HERP) 

Most  Navy  platforms  possess  systems  and  equipment 
which  radiate  energy  in  some  portion  of  the  electromagnetic 
spectrum.  The  radiation  can  constitute  a  hazard  for 
personnel .  The  degree  to  which  personnel  are  exposed  to  the 
radiation  depends  on  the  type  of  platform,  the  systems 
installed,  the  location  of  the  systems  (particularly  their 
antennas) ,  the  possibility  of  directive  antennas  illumin- 
ating areas  occupied  by  personnel,  etc. 

The  effects  of  exposure  to  this  radiation  can  be 
classified  into  two  categories,  direct  and  indirect.  Direct 
radiation  is  defined  as  the  energy  impinging  directly  on  the 
body.  Indirect  radiation  is  defined  as  the  capture  of  EM 
energy  by  a  metallic  object.  Resulting  voltage  and  currents 
cause  burns  when  the  metallic  objects  are  contacted  by 
personnel . 
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with  certain  reservations,  physiologists  generally 
agree  that  damage  to  personnel  from  direct  radiation  is 
primarily  a  heating  effect. 

In  recent  years,  however,  many  distinctly 
nonthermal  effects  have  been  documented,  some  of  which  have 
been  show  to  be  dependent  on  peak  powers  whose  average  value 
is  not  great  enough  to  produce  heating.  Frequency 
dependence,  with  no  heating,  has  also  characterized  many  of 
the  observed  effects.  While  the  full  significance  of  these 
effects  as  human  hazards  has  not  been  established,  the  fact 
that  they  occur  at  average  power  levels  considered  to  be 
negligible  suggest  that,  at  the  least,  an  awareness  of  their 
existence  should  be  assumed.  Some  recorded  nonthermal 
effects  listed  in  MIL-HDBK-238 (NAVY)  1973  are: 

(1)  Minor  changes  in  human  blood  properties  upon  exposure 
to  EM  energy  of  proper   frequency  and   intensity. 

(2)  Auditory  response.  Certain  people  hear  a  buzz  when 
exposed  to  microwave  radiation.  The  sensation  of 
sound  is  probably  not  the  microwave  frequency  but 
response  to  the  pulse  repetition  frequency. 

(3)  Abnormalities  of  the  chromosome  structure  occurring 
upon  exposure. 

(4)  Movement,  orientation,  and  polarization  of  protein 
molecules  in  pulse  RF  fields. 


61 


(5)  Unexplained   response   of  man   to   radar.   Epigastric 

distress,    emotional   upsets,     and    nausea   may 

occasionally  occur  at  as   low  as  5-10  mW/cm^    and  are 

most   commonly   associated  with  the   frequency   range 

from 

8  X  10^  to  12  X  10^  MHz 

Changes  in  the  transport  rate  of  materials  across  the 

blood-brain  barrier  in  humans  and  animals  have  also  been 

observed.    The  significance  of  these  changes  is  not  yet 

understood. 

The  Navy  RADHAZ   standard  follows  the  current  ANSI 

standard  for  3  to  3  0  MHz  with  constant  values  below  3  MHz. 

Personnel  shall  not  be  exposed  to  a  power  density  which, 

when  averaged  over  any  0.1  hour  period,   exceeds 

900/ f  (MHz)  mW/cm'  in  the  frequency  domain  of  3  to  3  0  MHz. 

Neither  the  root  mean  squared  electric  field  strength  (E) 

nor  the  root  mean  squared  magnetic  field  strength  (H)  may 

exceed  the  following  values  when  averaged  over  any  0.1  hour 

period: 

E  =  -J^^ —  V/M 
f  MHz    ^/^ 

H  =  ^'11 —   A/M 
f  MHz     ^ 

(These  are  the  electric  and  magnetic  field  strengths  roughly 
corresponding  to  electromagnetic  waves  in  free  space  to 
which  the  value  of  power  density  show  above  may  be  as- 
signed.) For  a  condition  where  exposure  is  not  regular  in 
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time  or  continuous  in  level  over  the  0.1  hour  period,  the 
equivalent  energy  fluence  level  of  one-tenth  of  the  power 
density  number  in  units  of  mW-h/cm^  may  be  use  as  the  limit 
of  exposure  for  any  0.1  hour  period.  In  situations  where 
measurements  of  two  or  more  quantities  are  available,  the 
most  restrictive  shall  be  used  as  the  limiting  factor. 

RF  burns  constitute  an  indirect  radiation  effect. 
Ship  communication  transmitters  induce  voltages  on  various 
metallic  structures  such  as  underway  replenishment  gear, 
booms,  and  loading  hooks.  The  voltages  are  capable  of 
causing  painful  burns  to  personnel  who  come  in  contact  with 
these  structures. 

It  has  been  established  that  140  volts  or  greater, 

measured  from  the  object  to  ground,  can  cause  burns  to 

persons  touching  the  object. 

b.   Hazards  of  Electromagnetic  Radiation  to  Ordnance 
(HERO) 

Hazards  of  electromagnetic  radiation  to  ordnance 

(HERO)   stem   from  the  use  of  sensitive  electroexplosive 

devices  (EEDs)  in  ordnance  systems.   The  principal  emphasis 

is  on  protecting  electroexplosive  devices  such  as  squibs, 

detonators,   explosive  switches,   and  ejection  cartridges. 

Premature  actuation  of  some  of  these  devices  can,  of  course, 

have  dire  consequences,  such  as  rocket  ignition  or  warhead 

detonation  on  a  flight  deck.   Premature  actuation  by  RF  can 

also  impact  the  reliability  of  a  system,  as  for  example  when 

the  power  supply  for  a  missile  guidance  system  is  expanded 
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inadvertently  by  the  RF  initiation  of  a  squib  switch.   In 
addition,   EED  firing  characteristics  can  be  altered  by 
electromagnetic  induced  heating. 

For  the  above  reasons  NAVMAT  Instruction  5101.1 
requires  that  weapons  systems  and  devices  containing  EEDs  be 
reviewed  and  tested  (if  deemed  necessary)  and  positive 
certification  obtained  that  they  can  be  handled  with 
impunity  in  the  maximum  predicted  electromagnetic  environ- 
ment before  introduction  for  service  use.  Three  classifica- 
tions pertinent  to  HERO  for  ordnance  items  have  been 
established.  The  classifications  are  based  on  weapons 
susceptibility.  The  degree  of  susceptibility  is  dependent 
upon  the  electromagnetic  environment,  the  potential  for  the 
induction  of  electromagnetic  energy  into  the  ordnance 
system,  and  the  characteristic  of  the  EED.  Items  that  are 
negligibly  susceptible  and  require  no  field  intensity 
restrictions  beyond  general  requirements  during  all  phases 
of  normal  employment  are  classified  as  HERO  SAFE  ORDNANCE. 
Items  that  are  moderately  susceptible  and  require  field 
intensity  restrictions  for  at  least  some  phases  of  employ- 
ment are  classified  as  HERO  SUSCEPTIBLE  ORDNANCE  SYSTEMS. 
Items  that  are  highly  susceptible  and  require  severe  field 
intensity  restrictions  for  some  or  all  phases  of  employment 
are  classified  as  HERO  UNSAFE  ORDNANCE. 

An  extensive  testing  program  exists  to  define  HERO 
problem  areas  and  ensure  the  safety  and  reliability  of 
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ordnance  items.    Tests  are  conducted  in  the  maximum  RF 

environment  to  which  ordnance  items  will  be  exposed  during 

the  stockpile-to-launch  sequence.    This  implies  that  HERO 

SAFE  ORDNANCE  items  have  been  tested  to  be  safe  and  reliable 

in  the  criteria  environments. 

c.   Hazards  of  Electromagnetic  Radiation  to  Fuels 
(HERE) 

The  possibility  of  accidentally  igniting  gasoline 

at  shore  facilities  and  aboard  aircraft  carriers  and  other 

ships  handling  fuel  has  been  considered  when  rf-indu«  -^d 

sparks   have   been   observed.     Many  years   ago   the   Navy 

conducted  extensive  tests  and  investigations  of  ignition  of 

fuel  vapors  by  rf-induced  sparks  and  was  found  that  the  risk 

of  ignition  was  small  because  the  following  conditions  must 

occur  simultaneously  for  ignition  to  take  place: 

(1)  Fuels  must  be  heated  above  their  flash  points. 
Fuels  will  not  ignite  unless  oxygen  is  provided  in 
certain  specific  and  exact  proportions.  Fuels  will 
also  not  detonate  unless  certain  exact  proportions  of 
oxygen  and  fuel  are  maintained. 

(2)  Sufficient  energy  must  be  provided  to  sustain 
ignition.  Although  radiant  energy  is  the  primary 
source  of  energy  associated  with  HERF  phenomena,  no 
mechanism  exists  by  which  this  energy  can  interact 
directly  with  the  fuel.  Radiated  energies  are  capable 
of  causing  arcing  in  susceptive  structures.  It  is 
these  arcs  which  produce  ignition  of  the  fuels. 
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An  arc  is  a  voltage  breakdown  between  two  fixed 
electrodes.  It  takes  a  certain  amount  of  voltage  to  break 
down  the  dielectric  between  the  electrodes,  but  it  takes 
energy  behind  this  source  of  voltage  to  sustain  the  arc  for 
sufficient  duration  to  ignite  fuel  mixtures.  A  spark,  on 
the  other  hand,  is  formed  when  two  metallic  objects  in 
electrical  contact  are  separated  and  there  is  sufficient 
energy  between  these  electrodes  to  sustain  the  spark  and 
thus  cause  the  ignition.  This  is  similar  to  the  situation 
in  which  a  fueling  nozzle  is  withdrawn  from  the  tank  opening 
when  an  aircraft  is  being  fueled.  The  radiation  field 
intensity  in  the  proximity  of  a  radiating  antenna  could  be 
sufficient  to  induce  the  required  energy. 

From  actual  measurements,  it  has  been  determined 
that  a  spark  of  energy  of  50  V-A  is  required  to  ignite 
gasoline  in  an  explosive  vapor  test  device.  Recently,  some 
attempt  has  been  made  to  relate  the  fuel  hazard  to  electric 
field  intensity,  the  primary  result  has  been  to  show  that 
the  igniting  electric  field  intensity  is  a  function  of 
frequency  and  is  minimal  in  the  upper  hf  band. 

In  recent  years,  the  probability  of  occurrence  of 
this  particular  hazard  has  been  further  decreased  with  the 
following  advances  of  technology: 

(1)  Ship  Design.  In  the  design  of  ship  topside 
arrangements,  the  location  of  fuel  handling  stations, 
fuel   vents,   etc,   is   considered   relative   to   the 
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placement  of  RF  transmitting  antennas.  Every  effort 
is  made  to  locate  these  antennas  away  from  the  fuel 
stations  and  vents. 

(2)  Fuel  Nozzle  Modifications.  Heat  shrinkable  tubing  is 
applied  to  nozzles  used  in  conventional  gasoline 
fueling.  This  prevents  metal-to-metal  contact  between 
the  nozzle  and  the  filler  tube,  thus  eliminating  the 
possibility  that  a  spark  will  occur.  Another 
hazard-reducing  development  in  the  aircraft  refueling 
area  has  been  the  introduction  of  the  pressurized 
fueling  approach.  In  this  procedure,  the  fuel  nozzle 
contacts  the  metallic  filler  tube  of  the  aircraft 
prior  to  the  opening  of  the  valve  to  the  tank.  Thus, 
if  a  spark  should  occur,  there  would  be  no  flammable 
fuel-air  mixture  to  ignite. 

(3)  Fuel.  The  type  of  fuel  now  being  most  extensively 
used  is  JP-5,  which  has  low  volatility  compared  to 
gasoline.  Because  JP-5  does  not  have  the  volatility 
of  gasoline,  it  is  not  likely  to  produce  flammable 
fuel-air  mixtures  under  accidental  circumstances,  this 
has  done  much  to  reduce  the  fuel  hazard  aboard  ships. 
The  fuel  oil  which  powers  the  ships  is  still  less 
flammable  than  JP-5. [Ref .  4: pp.  4-9] 
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